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The use of complex analysis for computing one-loop scattering amplitudes is naturally induced by generalised 
unitarity-cut conditions, fulfilled by complex values of the loop variable. We report on two techniques: the cut- 
integration with spinor-variables as contour integrals of rational functions; and the use of the Discrete Fourier 
Transform to optimize the reduction of tensor-integrals to master scalar integrals. 



1. GENERALISED UNITARITY 

The application of unitarity as an on-shell 
method of calculation [l|2j is based on the prin- 
ciples that products of on-shell tree-level ampli- 
tudes produce functions with the correct branch 
cuts in all channels; and that any one-loop ampli- 
tude is expected to be expressed, by Passarino- 
Veltman reduction, as a linear combination of 
scalar master integrals (MPs), that are charac- 
terised by their own, leading and subleading, sin- 
gularities [3]. 

Dimensionally-regulated amplitudes can be de- 
composed in terms of MPs in shifted-dimensions 
as, 



A 



(4-2e) _ 



E 

n, 0<j<j„ 



c nj x M^t 2t+2]) 



(1) 



where A 



(4-2e) 
JV 



is any ./V-point amplitude in D- 
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dimensions (being D 
is a n-point MPs in (4 — 2e + 2j)-dimensions, with 
n e {0, 2,3,4,..., N}, and j max depending on the 
process. In Eq.|T]), the coefficients c n j do not de- 
pend on D, and the whole Independence is em- 
bedded in the definition of M 's |4|2j , 
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where J\. is the generalised solid-angle, and 
In 2t+2 ^ is the scalar n-point function in shifted 
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dimensions. The decomposition (fTJ) could be fur- 
ther simplified with the help of recurrence rela- 
tions linking higher-point integrals to lower-point 
ones [5]. From Eq.lJT]), it is clear that the com- 
putation of the amplitude requires the knowledge 
of two types of ingredients: the MPs, and their 
coefficients. In the following, we focus on the de- 
termination of the latter ones. 

The principle of a unitarity-based method is 
the extraction of the rational coefficients, c nj -, 
by matching the multiparticle cuts of the ampli- 
tude onto the corresponding cuts of the MPs. By 
considering 6 2 7J the splitting of Ld, the loop 
momentum in (4 — 2e)-dimension, into its four- 
dimensional component, L, and its orthogonal 
complement i(_ 2e ), as Ld = L + L(_ 2e ) (with 
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dimensional integration measure can be written 
as a convolution of a four-dimensional integration 
and an integration on /i 2 , 
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By applying the above splitting to both sides of 
Eq.{T]), the four-dimensional kernel of the ampli- 
tude, A™ , can be read as expressed in terms of 



(4) 



four-dimensional n-point MPs, I, 
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We notice that the /^-dependence of A^ is due 
to the presence of [x 2 in all the denominators of 



as additional mass-term, and to the poly- 



tion on fj, , 



nomial coefBcients, c n j x (/i 2 )-'. It is therefore 
evident that to find out c n j , it is sufficient to com- 
pute the /^-polynomials [14] that are the coeffi- 
cients of the four-dimensional Mi's, In \ in the 
decomposition ([6]). 

1^ are functions determined by their own 
branch-cuts. Generalised unitarity is a very ef- 
fective tool to extract the rational coefficients 
of functions by exploiting their singularity struc- 
ture, which is accessed by imposing (on-shellness) 
cut-condition to propagating particles, 

{q 2 - m 2 + iO)- 1 -► (2ttz) <5 (+ V - m 2 ) . (7) 

In general, the fulfillment of multiple-cut condi- 
tions requires loop momenta with complex com- 
ponents. Since the loop momentum, L, in 
Eq.© has four components, the effect of the cut- 
conditions is to freeze some of its components, 
when not all, according to the number of the cuts. 
With the quadruple-exit [5] the loop momentum 
is completly frozen, yielding the algebraic deter- 
mination of the coefficients of J„ ,(n > 4); the 
spinorial integration of the double-cut [91101111375] 
and triple-cut [15116117) lead to the reconstruc- 



tion of 1% - and /^-coefficients; while the co- 
efficients of /q 4 ' are detected by single-cut. In 
cases where fewer than four denominators are 
cut, the loop momentum is not frozen: the free- 
components are left over as integration variables 
of the phase-space. We will discuss two strate- 
gies for dealing with the degrees of freedom rep- 
resented by those variables: i) analytic inte- 
gration of the phase-space with spinor variables 
[9110111112113114] : ii) algebraic decomposition of 
the integrand by means of the Discrete Fourier 
Transform (DFT) [23]. 
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2. DOUBLE-CUT AND SPINORS 

By using the splitting of the loop variables as 
above, the D-dimensional double-cut in the P- 
channel of any amplitude, i.e. the double-cut of 
the l.h.s of Eq.fTJ), can be written as a convolution 
of a four-dimensional double-cut and an integra- 
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is the double-cut of the l.h.s of 
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Eq.©, 

A(A<$) = J d^L S^(L 2 -M 2 -^ 2 ) 

x5 i+ \{L-P) 2 -M 2 -p 2 ) A^ cc A| ree , (9) 

with A\ ree being the tree-amplitudes sewn in the 
cut. We found it convenient to decompose [11] 
the four-dimensional loop variable, L, in terms 
of a massless momentum £, and the momentum 
accross the cut, P, 

L v = t £ v + z Q P v , (10) 
t=(l-2z )P 2 /(£\P\£] , (11) 
z = {P 2 + M 2 - M 2 - y/\ - 4/i 2 )/2P 2 , (12) 
A = (P 2 ) 2 + (Af 2 ) 2 + (M 2 ) 2 + 

-2P 2 M\-2P 2 Ml~2MlMl , (13) 

with zq being the anomalous threshold, and A, 
the Kallen function. With the above transforma- 
tion, A(A N ) can be written in terms of spinor- 
variables [9|10j . \£) and \l] (associated to the 
massless momenta, £, through (. = and can 

be cast as a sum of terms whose general structure 
reads, 

A(^ 3 ) = ]T A, , Ai = hi d£)[£ d£] li ,(14) 
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where Pi and Pi can either be equal to the cut- 
momentum P, or be a linear combination of ex- 
ternal vectors; and where the pi's depend solely 
on one spinor flavour, say \£) (and not on \£]), and 
may contain poles in \£). We give as understood 
the dependence of Ii on p, 2 , through the variable 
zq. The explicit form of the vectors Pi and P 2 
in (|15[) is determining the nature of the double- 
cut, logarithmic or not, and correspondingly the 
topology of the diagram which is associated to: if 
Pi = P2 = P, 



li= P i(\£))[vl} n /(t\P\t] 
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and the result will be non-logarithmic, hence cor- 
responding to the cut of a 2-point function with 
external momentum P; if Pi = P, P-2 ^ P or 
Pi =/= P2 =/= P, one proceeds by introducing a 
Feynman parameter, to write li as, 

v„ pi(\e})tn£\ n 



li = (n + 1) / 
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(e\R\e] n + 2 

fL = zf> 1 + (l-x)j>2, (18) 
and (because of the parametric integral) the re- 
sult is logarithmic, hence containing the cut of 
a linear combination of n-point functions with 
n > 3. The spinorial structure of Eq. (flU]) and 
Eq. (fl7|) is the same. Therefore, we discuss the 
spinor integration of the latter, because it is more 
general. 

2.1. Contour Integrals 

One can proceed with a change of variables [TH] , 
decomposing \£) and \£] into two arbitrary mass- 
less momenta, p and q (light-cone decomposition) , 

Vp,q:q 2 =p 2 = , (19) 

\t) = \p)+z\q) \e] = \ p ] + z\ q ] (20) 

(£ dt)[l dl] = -{q\p\q] dz dz . (21) 

Its effect on A* reads, 
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X n+2 (z,z) 
X(z, z) = (p\R\p] + z(q\R\p] + 

+z(p\R\q] + zz(q\R\q] . (23) 

One may observe that the z-2-integrand can be 
written as a total derivative with respect to z 
r-i 
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with ^(2) = ((p\P\q] + z(q\P\q\), so that the spinor 
integration has been turned into a contour inte- 
gral of a rational function in z, 

A, - / dx (1 - x) n (q\p\q\ 
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The z-integral can be performed by Cauchy's 
residue theorem, summing the residues at the 
poles in z (substituting as well z — z*). There 
are two sources of poles to account for: i) the 
poles contained in Pi(z); ii) the pole due to £(z), 
whose value is 
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To complete the integration of A; in (|25j) . one 
has to perform the parametric integration which 
is finally responsible for the appearence of loga- 
rithmic terms in the double-cut. On the contrary, 
the spinorial integration in (fT6|) would generate 
a contribution without branch-cuts. We remark 
that the role of \£) and \£] in the integration could 
be interchanged. 

At the end of the phase-space integration, by 
adding up all the A^'s, one finally gets a result 
whose structure is 



a(A 4) )= E 



x (P 



2 v 



(27) 



2<n,j 



corresponding to the double-cut of Eq.©. Out 
of ([27]) . it is possible to extract the polynomial 
coefficients, c n j x (p 2 y (n > 2): the coefficient 
of Iq 4 \ the tadpoles, cannot be detected within 
the double-cut, and their determination should be 
provided by independent informations on the am- 
plitude. We recall that the /^-dependence of the 
coefficients originates from the understood pres- 
ence of zq, given in Eq.JTJ]). 

We observe that a proper choice of the mo- 
menta p and q, entering the change of variables 
(|20jl . can simplify dramatically the calculation. 
For instance, they determine the value of the 
z^-polc, given in Eq. (|26[) : given q^, and the 
cut-momentum P M , the choice p^ = P^ — q^ x 
P 2 /(q\P\q], would yield z^ = 0. 

The phase-space integration just discussed was 
used succesfully for an analytic computation of 
non-trivial one-loop corrections. In particular, its 
four-dimensional (massless) version [9 10 16J has 
been applied to complete the non-supersymmetric 
cut-constructible term of the six-gluon amplitude 
in QCD 10], to compute the six-photon ampli- 
tude in QED [20119] . and the cut-constructible 
term of a general MHV amplitudes involving a 



Higgs plus 7i-gluons in QCD (in the heavy-top 
limit) [21]. 

Recently, the efficiency of spinor integration 
has been pushed to achieve closed analytic forms 
for the generating formulas of the coefficients of 
(2 < n < 4), for an arbitrary massive pro- 



cess [13|14] , which together with 1^ constitute a 
basis of functions in four-dimensions, hence, due 
to the relation among Eq.© and Eq.lJTJ), in D- 
dimensions. The formulas presented in [13ll4j - 
too long to be shown here - can be evaluated, 
without performing any integration, by specializ- 
ing the value of input variables that are specific to 
the initial cut-integrand as assembled from tree- 
level amplitudes. 

We have as well recently released the pack- 
age S@M (Spinors @ Mathematica) [22] , that imple- 
ments the spinor-helicity formalism in Mathemat- 
ica. The package allows the use of complex-spinor 
algebra along with the multi-purpose features of 
Mathematica, and it is suitable for the algebraic 
manipulation and integration of products of tree 
amplitudes with complex spinors sewn in gener- 
alised unitarity-cuts. 

3. OPTIMIZED REDUCTION 

As an alternative to any phase-space integra- 
tion, in 24 25j there was proposed a very efficient 
method for the reconstruction of the coefficients 
in the decomposition ©. In what follows, I limit 
the discussion to the so called cut-constructible 
term of a scattering amplitude, that corresponds 
to the poly-logarithmic structure arising when 
Eq.© is evaluated at ^l 2 = 0. I will sketch 
the reconstruction of the complete /^-dependence 
[26127128129] at the end of the section. The by- 
now known as OPP-reduction allows the numerical 
reconstruction of c„o, by solving a system of al- 
gebraic equations that are obtained by: i) the 
numerical evaluation of the integrand at explicit 
values of the loop- variable, on the one side; ii) 
and the knowledge of the most general polyno- 
mial structure of the integrand itself [30j , on the 
other one. The values of the loop momentum used 
for the numerical evaluation of the integrand are 
chosen among the set of solutions of the multiple- 
cut conditions, i.e. the solutions of the system of 



equations obtained by imposing the vanishing of 
the cut-denominators. 

3.1. DPP-Reduction 

The starting point of the OPP reduction 
method 24 25] is the general expression for the 
integrand of a generic m-point one- loop amplitude 
that can be written as 



A m (q) = 



N(q) 
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where N(q) is the four-dimensional numerator of 
the amplitude]! The main formula of the 0PP- 
algorithm is the expression of N(q) in terms of 
the denominators £),-, 



N(q) = £>„(<?) 
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with 
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By inserting (|3"0|) back in l[28p. one exposes the 
multi-pole nature of A rn . The coefficients of 
the multi-pole expansion can be further split 
in two pieces: a piece that still depends on q, 
parametrized by d, c, 6, a, that vanishes upon in- 
tegration, and a piece that does not depend on 
q, parametrized as d, c, b, a. Such a separation is 
always possible, as shown in [24] . and, with this 
choice, the latter set of coefficients corresponds 
to the ensemble of the coefficients of In(fi 2 = 
0) , (ne {0,2,3,4}): a, b, c, din ([30]) correspond 
respectively to c o, c 2 o, c 30 , C40 in (BJ). 



A m (q) is the integrand of (fi? = 0), denned in Eq.JBJl 



3.2. Top-Down System 

The goal of the algorithm is reduced to the alge- 
braical problem of fitting the coefficients d, c, b, a 
by evaluating the function N(q) a sufficient num- 
ber of times, at different values of q, and then 
inverting the system. Accordingly, let us define 
the following functions, 
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R' ijh {q) ee (N(q) - A 4 (g))( ]J Dp) 1 , (36) 

4 m — 1 

R'r(q) ee (N(q) - ]T A a (q))( JJ Dp)' 1 , (37) 

4 m — 1 

<(g) = (JV(g) - ]T A„(«))( I] • ( 38 ) 

We as well define as {q}ij k e the set of the solutions 
of Di = Dj = D k = D( = 0. Having defined our 
setup, from Eq. ([30|) we can derive the following 
sets of equations: 



Rijke(q) = dijki + dijke(q) 
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(39) 
(40) 
(41) 
(42) 



which must be solved necessarily in cascade, top- 
down: in Eq. ((39|) . N(q) is a known quantity, 
namely an input of the algorithm; but the l.h.s 
of each other equation becomes a known quan- 
tity (numerically evaluable), only after solving 
the equation which preceeds it. 

3.3. Polynomial Structures and DFT 

An important observation is due. The r.h.s of 
each of the equations (|39p - (l42[) is a polynomial 
function. Without presenting their explicit ex- 
pressions (see 23 for the detailed presentation), 
the general structure is the following: the vari- 
ables are the components of q not-frozen by the 
cut-conditions; the degree is known; while the 



coefficients are the unknowns to be determined. 
The problem to be tackled is thus a well known 
mathematical subject: polynomial interpolation. 
In order to find out the coefficients of a polyno- 
mial, one can avoid the numerical inversion of a 
system, which is a very delicate operation, due 
to the possibility of a vanishing determinant in 
critical kinematic regions. 

The Discrete Fourier Transform (DFT) is a 
very efficient tool to extract the coefficients of 
a polynomial, by evaluating it at special values 
of the variables [13 23 31 . Let us show how it 
works in the case of a polynomial of degree n in 
one variable, x, defined as, 



(43) 



At the first step, one generates the set of discrete 
values P n , k [k = 0, n), 



Pn,k = Pn{xk) Cg P l e 



(44) 



by sampling P n {x) at (n + 1) equidistant points 
on the p-circle, 
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(45) 



At the second step, using the orthogonality 

n-l 

e 2m i i e~ 2 ^ i = n 6 kk , , (46) 

3=0 

one can obtain the coefficient eg simply by pro- 
jection, 



n + 1 f-^ 



(47) 



k=0 



In fact, the r.h.s of Eq. (l42l) is a degree- 1 poly- 
nomial in a single variable, whose coefficients are 
easily determined by the semi-sum and the semi- 
difference of two numerical values of R. But the 
r.h.s of Eqs. (|40j) - (j42|) are multivariate polynomi- 
als of higher degree. To find out their coefficients 
we used a modified DFT, that is a Fast Fourier 
Transform-like algorithm, suitable to minimize 
the number of the numerical calls respectively of 
R',R", and R'", being exactly the same as the 



number of the unknowns, and to avoid the kinc- 
matical singularities emerging at the vanishing of 
the circle-radius p. In so doing, one can deter- 
mine all the unknown coefficients, among which 
the O th -order ones, respectively dijke, Cijk, hj, a>i, 
correspond to the coefficients of the Mi's in four- 
dimension. 

For the reconstruction of the complete [in- 
dependence of the coefficients in Eq.©, the de- 
composition (|3"0)) . must be slightly extended to 
account for the presence of p, 2 [28129] . The start- 
ing point, in this case, is A${p?), which con- 
tains a numerator N(q, p 2 ) and denominators 
L>i = Di — p 2 . The reduction proceeds exactly 
as above, with the difference Eqs. ([55)1 - (|4"2"1) con- 
taining an extra dependence on p 2 . Since the 
/^-dependence is still polynomial, one can use 
the DFT also in this case, having to deal with 
R, R' , R", and R" 1 with p 2 as additional variable 
[13] . The flexibility of the projection-procedure 
hereby presented extends its range of applicability 
to tackle the determination of the coefficients of 
polynomial structures wherever should this issue 
occur. We finally remark that the parametriza- 
tion of the free (integration) variables as com- 
plex unitary phases yields as well a very effective 
performance of Cauchy's residue theorem within 
the contexts of factorization- and unitarity-based 
methods, where the on-shellness properties are 
naturally captured by polar structures of complex 
phases. 
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